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In a resent paperS Eriksen and Hedegard (EH) claimed that they "clarified" our path integral 
calculationcl and found an error in our analysis. Here we will demonstrate that this claim is, based 
on mathematically incorrect and physically meaningless procedure invented by the authored. The 
main problem in the EH analysis is due to their erroneous calculation of Gaussian path integrals 
and has nothing to do with the momentum recoil of a particle interacting with an external bath. 



In a very recent preprint!! Eriksen and Hedegard (EH) attempted to reanalyze a part of our calculation! and 
"demonstrated" that their version of the effective action for an, electron interacting with the rest of the fluctuating 
' electronic environment is different from one found in our papern. This difference - according to EH - should explain 
disagreement of our results with ones found by Aleiner, AJtshuler and Gershenzon (AAG)El and invalidate our analysis. 

On one hand, we welcome an attempt of the authorsa to go deeper into our path integral analysis.. As a result of 
that our work was upgraded from "profoundly incorrect calculation" El to beinc' "not without merits"EI. In some sense 
■ EH also defended our work indicating that the previous "vague statements"! 3 ] concerning "problems in semiclassical 
calculations when classical orbits, intepsect" have little to do with our calculation. We fully agree with that. 
t-H , On the other hand, the resultso andta apparently differ at low temperatures and according to both AAG and EH the 
problem should be found in jour calculations. This is an erroneous opinion. Inu we will demonstrate that the results 
of our path integral analysisa axejully equivalent to those of AAG at the level of the Golden-rule-type perturbative 
analysis employed by the authorsa. In other words, if we expand our path integral in powers of the effective action, 
keeping only the first term of this expansion we arrive exactly at the AAG resultsQ. The actuaLdifference betweentl 
andtJ is that we do nq± expand our path integral and go further than the Golden rule approach! 3 ]. This issue will be 
discussed in details inta. n 
ON ■ Here, however, we are concerned with a purely technical problem. Namely we will analyze the clainu concerning 
the error in our path integral aaalysisa. EH proceeded quite far in checking our calculations and arrived at our Eqs. 
(43)0 (equivalently Eqs. (11) intl) where - according to them - the problem occurs. Although EH did not point out 
{""j where a mistake can appear in our analysis they presented their own derivation which - according to them - should 
1 , \ "repair" our calculation. Already at this point we would like to note that the part of our derivation which EH would 
like to repair is the exact quantum mechanical calculation involving no approximations. It yields the effective action 
which - as it was already pointed out - on a perturbative level allows to obtain the results identical to those of AAGcl. 
O . But it is namely this effective action which is challenged by EH. _ 

Since we cannot exclude that EH simply misunderstood the corresponding part of our derivationB we will reproduce 
• i-H , it again carefully explaining some of our steps in order to avoid possible further misunderstandings. Then we analyze 
1 the derivation suggested by EH. 

The first point where EH do not agree with our results is the path integral reperesentation of the electron evolution 
■ - - ■ operators 
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-i l dt'H ia {t') 
t 



(1) 



where 



Hi( P , r) = 7^ + U(r) - eV+(t, r) - i [l - 2n(H (p, r))]eV~(t, r), 

H 2 (P, r) = |- + U(r) - eV+(t, r) + \eV~(t, r) [l - 2n(H (p, r))] , (2) 

EH correctly pointed out that the operators n(Ho(p, r)) and V~(t, r) do not commute and this makes them different 
from the ordinary single particle Hamiltonians. However, the path integral representation of the evolution operators 
(ffl) always has the following standard form (see e.g.El): 
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Ui, 2 {t 1 ,t 2 ;r f ,r i ) = J Vr(t') J Vp(t') exp i J dt' (pr - Hi, 2 (p,r)) 
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The commutation rules need to be specified separately. E.g. considering the expression for the operator U\, one 
can introduce an infinitesimal time shift (+0) and to write the product (1 — 2n)V~ in the action as [1 — 2n{p(t + 
0), r(t))]V~ (r(t)). This infinitesimal time shift guarantees that in the perturbation expansion the operator 1 — 2p 
will be always placed before the operator V~ . The expression (|^) together with the commutation rules is the exact 
quantum mechanical formula. 

The next step of our analysis is to calculate the kernel of the electron evolution operator J in the presence of 
interaction with the effective electromagnetic environment. This function can be expressed in terms of the path 
integral 

ri(t)=r lf r 2 (t)=r 2f 

J(t,t'- 1 r 1 f : r 2 f,r 1 i : r2i)= J Vri J Vr 2 J VpiVp 2 x 

ri(t')=Tn r 2 {t')=r 2i 

p iS , o[T- 1 ,p 1 ]-iS , o[r 2 ,P2]+'i /* dt" f dr(f-V++f+V-) 
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where the electron action Sq[x,p] has the form 

S [r,p] = J dt"(pr-^-U(r)\ 

f 

and the "charge densities" f~, f + arc defined by the equations: 
f-(t,r) = eS(r- ri (t))-e5(r-r 2 (t)), 

/+(*, r ) = Ue[l- 2n( Pl (t + 0), n(t))] S(r - n(t)) + e[l - 2n(p 2 (t + 0),r 2 (t))]S(r - r 2 (t))^j . 
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Averaging over the fluctuating fields V + , V in Eq. (|4|) amounts to calculating Gaussian path integrals with the 
action 

dud 3 k M fc 2 e(uvfc) T/+ 

(2tt)4 
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f dujd 3 k Tr . , , fc 2 Ime(w, fc) , / oj \ Tr . , . 
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where e(w, k) is the dielectric susceptibility of the system. This averaging can be easily (and, of course, exactly!) 
performed and yields 

ri(t)=r lf r 2 (t)=r 2} 



J(t,t';r 1 f,r 2 f,rii,r 2i )= J Vri J Vr 2 J VpiVp 2 x 



r\{t')—r\i r 2 (t t )—r 2 i 

x exp {iS [ri,pi] -iS [r 2 ,p 2 ] - iS R [r 1 ,p ll r 2 ,p 2 ] - Si[n, r 2 ]}\ (8) 

where 
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SAr Xl r 2 \ = e — J dti J dt 2 (l(t 1 -t 2 ,r 1 (t 1 )-r 1 {h))+I(t 1 -t 2 ,r2{ti)-r 2 (t 2 ))- 
t> v 

- * 2> r 1 (ti) - r 2 (t 2 )) - I{h ~ t2,r 2 (h) - nfe)) J. (10) 
The functions i? and / are defined by the equations 

/ffK^M£>""* (12) 

The combination in the exponent (||) is the effective action. It is defined on the Keldysh contour, the last two 
terms Sr and Si depend on two coordinates r\ and r 2 and are defined by nonlocal in time expressions |Io|) . The 
influence of the environment, the Pauli principle and the commutation rules between p and V~ are fully accounted 
for. Proceeding from ([!]) to ( |12] ) we made no approximations. The quasiclassical approximation is used already at a 
later stage of our calculation. Then the infinitesimal time shift does not play any role and can be ignored. But since 
EH challenge our effective action, namely its part (||), we will not discuss further steps of our calculation^. 

Let us now examine the derivation suggested by EH. They define the quantity 

1 = (Pf, x f\Ui(t)\x l ,p l ) v ± (13) 

(their Eq. (20)) which - although can be considered as a mathematical object - has little physical meaning in our 
problem. Note that U\ here is the evolution operator defined only on the forward part of the Keldysh contour, while 
both the fields V and the effective action are (and should be!) defined on both forward and backward part of this 
contour. Already here EH have a problem: how is it possible to derive the effective action in the presence of the 
environment working only with one part of ithe Keldysh contour? This contradics to all what is known about this 
matter. But this is a textbook issue (see e.g.H) and we are not going to discuss it in more detais here. 

Since EH would like to calculate the quantity I (|l^) we will first demonstrate how to do it correctly. If one uses 
the exact quantum path integral representation (|l|) of the operator U\ (t) , its kernel can be defined as 

r(t)=x f 

J 1 (t,0;x f ,x i )= J VrJv P e tSo[r < p]+l Jo dt J dr ^ V++f " V ~), (14) 

r(0)=Xi 

with 

fr(t,r) = eS(r-r(t)), 

/+(*, r) = (e/2) [1 - 2n(p(t + 0),r(t))]*(r - r(t)). (15) 

The equations (|lj) and ( |l5|) should now be combined with (|l3|). The Gaussian path integral over the fields V + and 
V~ can be evaluated exactly and one finds 

r{t)—Xf 

1= J Vr J Vpe iSo ^ +iS "f^ r ' p \ (16) 
r(0)=sBj 

where the action S e f / is given by 

t t' 

S eff = e 2 J dt' J dt" (i(V + (t',r(t'))V + (t",r(t"))) 
o o 

+±[l-2n(p(t" + a)Mn)KV + (t'Mt'))V-(t'\r(t'%y (17) 



3 



This action depends on both coordinate r and momentum p. It does not contain any recoil! We repeat: the quantities 
/ (|l3| ) and S e f / (|l7j ) by itself have little physical meaning. But formally they can be defined and calculated as it was 
shown above. 

What do EH suggest instead of this simple exact calculation? 

To make the comparison with EH results easier, it is useful to rewrite the action fflTi) in the following form: 





Seff[r,p} = ^ dt dr + f+V~) (18) 



In order to evaluate their quantity / (jTj) EH used the procedure which is nothing else but the expansion of the 
exact expression (|l4| ) in powers of V r± up to the second order, averaging this first order correction over the fluctuating 
potentials and writing it again in the form of the exponent. Namely, their procedure reduces to the following set of 
approximations 



j _ { e iS [r,p]+i f*dt J dr(f-V + + f+V') 

i-i(/Wv- + /^-)) 



giSo [r,p] 



V^ 1 ,p,r 

t 



» (19) 

V^ 1 ,p,r 




Then they calculate the first order correction given by the last term in the square brackets and call the following 
construction "effective action": 

j Vp e iS [r, P ] /fjdtjdr {fiV+ + f+V- 

Seff M = 1 ^-——p— ^j—, — i 2( ) ! 

2 J Dp e lSo [ r 'P] 

This "action" is very different from the exact result ( |l7j , |l8| ) and is just slightly modified first order of the perturbation 
theory in interaction. EH do not even bother calculating this "action" and leave it in the form equivalent to fl2C| ) 
(see their eq. (27)). They evaluate it only on a straight line trajectory (their eqs. (28,29)). Apart from the fact that 
these trajectories are meaningless in the case of disordered metals at all relevant distances (i.e. distances longer than 
the mean free path), we do not understand how one can compare the action as a functional of all trajectories with a 
number calculated for one trajectory only. 

Anyway, in contrast to the EH claim the problem has nothing to do with the commutation relations but rather 
emerges from a wrong calculation suggested by EH. In order to illustrate that one can consider e.g. the muon motion 
which is not restricted by the Pauli principle. In this case the terms containing n(Ho(p,r)) in the Hamiltonians (Q) 
are absent and the problems with commutation of the operators do not appear at all. However, even in that case the 
exact action (|j!ljl8) remains very different from the action (p0|). This example demonstrates that the procedure used 



by the authorsu fails even in the simplest case. 

Now let us discuss the momentum recoil. EH concluded that our theory does not describe the effect of recoil because 
it is not observed in the action. This is wrong. Our expression for the effective action was obtained by the exact 
(within RPA) integration over all the degrees of freedom of the effective environment and hence includes all processes 
in all orders. It would be therefore very naive to try to interpret the full action only in terms of the simplest first 
order scattering processes. However, if one wants to describe only such processes one can easily do that proceeding 
perturbatively and expanding the path integral in powers of the exact action-, (§-^0|) . All correct recoils and other 
effects are fully reproduced. It was demostrated e.g. in Sec. 5 of our papera. There EH can easily find e.g. the 
combination n p ^k entering eqs. (90), (91), (93-95) or jfthe standard result from the standard many body theory" 
containing the combination "coth-tanh" (see eq. (97) ofl). But all that is only the simple Golden rule perturbation 
theory which is actually known to fail at low T. It is really hard to understand how the exact procedure can be 
"repaired" by making a rough approximation! 

One can assume that if one proceeds perturbatively in the interaction the expansion done by EH is legitimate. 
This is not true: although interaction may be small, the time integrals are large, and the corresponding terms should 
be always kept in the exponent in the long time limit. If one wants to expand one should keep all orders of this 
expansion. But to be on a safe side it is always better to calculate exactly as long as one can do it. This will never 
make things worse. 

Just for an illustration let us make a trivial exercise. Consider a Gaussian integral 
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Ii(a,b) — / dxexp(— ax 2 — 2bx) 



(21) 



In order to find out what is the leading dependence of I\ on b it is not even necessary to calculate this trivial integral. 
One just makes a shift and obtains I\ oc exp(6 2 /a). Then one can decide if one wants to expand the exponent or not. 
However if one expands in b first the result may be wrong depending on the value b 2 /a. This seems to be obvious. 
Consider now a slightly more complicated Gaussian integral 

N 

I N [a, b}= / dx n exp(-a„a; 2 - 2b n x n ) (22) 

^ n=l 

Analogously one obtains 

JV 

I N cx exp(J^ b 2 Ja n ) (23) 
n=l 

Again if the whole sum in the exponent is small one can expand. However, if one expands in all b n first and calculates 
the integrals over all x n one takes a risk of making a mistake. Some of 6^/a„ may be small, some others not and - 
especially if N is very large - the sum in the expoaent will not be small. 

EH's - and not only EH's - critique of our papercl is essentially on this level: "expand first, get-."coth-tanh" , obtain 
zero result at T — ► 0, compare with our resultH, observe the difference, conclude that our resultn is wrong". But we 
agree (and always knew it) that one gets "coth-tanh" if one expands first. What we are saying is that it is wrong to 
expand first. We welcome constructive critique of our results. But it is always better to understand what was actually 
done in the paper and criticize afterwards than vise versa. 

Summarizing, the EH claim that we made an error "when dealing with the real space path integral representation 
of quasiparticle propagation" is wrong. [By the way, we do not quite know what the quasiparticle in a disordered 
metal with interaction is, do EH know?] In contrast to EH beliefs our formalism does "allow-fbr tiny changes in the 
momentum and energy" . It actually does, it in all orders and not only in the first order like inlil. The erroneous result 
indeed occurs, it occurs in the EH papeiu and not because of "tiny changes" but rather because the authorsEl believe 
that it, is more accurate to calculate Gaussian integrals perturbatively than to do it exactly. The story with the EH 
papertl is over. 
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